Scattering matrices with block symmetries 



Karol Zyczkowski 

Instytut Fizyki M. Smoluchowskiego, Uniwersytet Jagiellonski, id. Reymonta 4, 30-059 Krakow, Poland 



00 



o : 
o ■ 

m ■ 

o : 
. 

o\ ■ 

>v 

i ' 

o : 
^ . 

o ■ 



13 



Scattering matrices with block symmetry, which corre- 
sponds to scattering process on cavities with geometrical sym- 
metry, are analyzed. The distribution of transmission coeffi- 
cient is computed for different number of channels in the case 
of a system with or without the time reversal invariance. An 
interpolating formula for the case of gradual time reversal 
symmetry breaking is proposed. 



I. INTRODUCTION 

Ensembles of random matrices introduced in context 
of the theory of nuclear spectra by Dyson H long time 
ago found a novel application in problems of chaotic scat- 
tering [^| D . The ^-matrix corresponding to time rever- 
sal invariant problems can be described by matrices of 
circular orthogonal ensemble (COE), while the circular 
unitary ensemble (CUE) is applicable if no antiunitary 
symmetry exists. 

The process of chaotic scattering in cavities with a ge- 
ometrical symmetry can be represented by an S'-matrix 
with block symmetry. Such matrices where recently in- 
troduced by Gopar et al. j|, who discussed both cases: 
with and without time-reversal symmetry. They found a 
link between the invariant measure of such ensembles and 
the measures of canonical circular ensembles and com- 
puted the distribution of transmission coefficient for one 
and two incoming channels. 

In this Brief Report we generalize the results of Gopar 
et al. || . Using an idea of composed ensembles of random 
matrices we simplify analytical calculation of the trans- 
mission coefficient T. This method allows us to obtain 
the distributions P(T) for interpolating ensembles corre- 
sponding to breaking of the time reversal symmetry and 
to treat the case of large number of channels. Generat- 
ing random unitary matrices of interpolating ensembles 
according to the method described in |Tc| l we verify nu- 
merically proposed distributions. 



II. BLOCK SYMMETRIC SCATTERING 
MATRICES 

We introduce the random S-matrices with block sym- 
metry following the method of Gopar et al. ||. The 
2M x 2M scattering matrix possess the structure 



S 



r t 
t r 



(1) 



where M x M matrices r and t describe reflection and 
transmission processes, respectively. S-matrix can be 
brought into a block diagonal form 



S' — RqSRq 



si 
s 2 



(2) 



where s\ = r + t and S2 = r — t are unitary matri- 
ces and i?o stands for rotation matrix consisting of M- 
dimensional unit matrices 1m 



R 



= J_fU 

V2 V- 1 



1m 
1m 



(3) 



Of great physical importance is the total transmission 
coefficient 



T = tr(ii t ), 



(4) 



which is proportional to the conductance G of the cavity, 
G — (2e 2 /h)T. Since t = (si — s 2 )/2, the transmission 
coefficient reads 



T 



M 
Y 



1 



Re[tr(si4)]- 



(5) 



Following Gopar et al. j^] we assume that unitary matri- 
ces si and S2 are statistically independent and pertain to 
the same universality class. For canonical circular ensem- 
bles of random matrices the joint probability distribution 
(JPD) of eigenphases is given by a single formula 



, <pm] = c p ,m TT I 



m 



(6) 



where Cp t M stands for normalization constant, while Up 
represents Poissonian, orthogonal and unitary circular 
ensemble for (3 equal to 0, 1 and 2, respectively. 

Formula (^) contains a product of two unitary matrices 
S1S2 (one can redefine the second matrix writing s' 2 = 
s\). Further calculation base on the relation linking the 
JPD of eigenvalues of the product of two unitary matrices 
drown independently from any canonical ensemble 



(.Up [<Pi,---, Pm] = p Up [pi-,---, <Pm] 



(7) 



The left hand side of the above equation formally rep- 
resents the JPD characteristic to spectra of a composed 
ensemble defined via product of two random matrices, 
each specified by a certain probability distribution. For 
unitary ensemble (/3 = 2) this equality follows from the 
invariance properties of CUE, which corresponds to the 
Haar measure on the unitary group. The same concerns 
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the case (3 = 0, since the circular Poissonian ensemble 
(CPE) can be defined by the Haar measure in the space 
of diagonal unitary matrices. In the case of symmet- 
ric unitary matrices ((3 = 1) this property follows from 
the definition of COE, which is invariant with respect to 
transformations s — ► s' = XsX T [fij, where X denotes 
any unitary matrix. Analyzed matrix S1S2 is similar to 

1/2 1/2 c . . , . . 1/2 , . , 

Si S2S1 ■ Since si is symmetric, so is , which may 
play the role of X in the invariance condition. The JPD 
of eigenvalues of S1S2 (averaged over the composed en- 
semble) is thus the same as this of of s 2 and is given by 
with (3 = 1. In spite of this result the measure of 
the composed ensemble containing products of two sym- 
metric unitary matrices differs form this characteristic of 
COE @. 

Since the JPD of eigenvalues does not define the entire 
probability distribution of an ensemble of random ma- 
trices, the formula (R) looses its meaning for non integer 
values of the parameter (3, characteristic to transition be- 
tween the canonical ensembles. However, in order to get 
possible interpolating formulae for distribution of trans- 
mission coefficients we will eventually allow (3 to take any 
real value in JO, 2]. 

Relation (0) allows one to write the transmission coef- 
ficient (j5[) in a simplified form 

T=^-iRe[tr(C/ /3 )], (8) 

and to obtain the distributions P(T) by averaging above 
formula over an appropriate ensemble of M x M unitary 
matrices Up. 

III. EXPECTATION VALUES 

Since (tr(Up)) = for any (3 one obtains 

M 

(T)p = y- (9) 

Let tr(Up) = z = ae vd . In a recent work |l3| the follow- 
ing average was derived for canonical circular ensembles: 
(a 2 )p = 2M/(2 + I3{M - 1)). Because the distribution 
of phases $ was shown to be uniform, the variances of 
both parts are equal: var(Re(z))=var(Im(z))=var(a)/2. 
Using (||) we get directly 

™(T h = { (T- { T h f) = 8 + ^ M _ iy (10) 

Observe that for (3 = 2 the variance equals to 1/8, in- 
dependently of the matrix size M, while for (3 = 1 one 
has var(T)i = M /A(M + 1), in accordance with earlier 
results [[J . The variance growths with a decreasing (3 and 
tends to M/8 in the limiting case (3 — *■ 0. 



IV. DISTRIBUTION OF TRANSMISSION 
COEFFICIENT P (T) 

A. The case M=l 

For M=l the "one dimensional matrix" U = e lVl and 
the phase ipi is uniformly distributed in [0, 2tt) for any 
ensemble. The variable t = cos^i) has thus the distribu- 
tion Pt(t) = l/[7T-v/ (1 — t 2 )]. The transmission coefficient 
equal in this case T = 1/2 — t/2 is therefore distributed 
according to 

for any value of (3. 

B. The case M=2 

Let us start deriving the distribution Pp{a) of the ab- 
solute value of trace a = \Tr(U)\. For M = 2 the JPD 
(0) reads 

P [ ( pi, i p 2 ]=Cp, 2 (sm<j>f, (12) 

where <fi = (cpi — ^2) /2. The module of the trace a = |1 + 
e i20| _ 2cos(/>, so employing ( p~2| ) we obtain the required 
distribution 

Pp(a) = c ^(4-a 2 )^, (13) 
where a € [0, 2] and the normalization constant reads 

In particular this distribution is flat for COE, while for 
CUE one gets a semicircle and recovers the result already 
mentioned in jl3| . According to (ph the transmission co- 
efficient equals T = 1 — (acosi?)/2. Due to rotational 
symmetry of the ensembles the distribution of phase i9 
is uniform in [0, 27r), so the distribution of t = cosd is 
P t (t) = l/[ny/(l-t 2 )]. Denoting x = at and using @ 
one writes an integral for the distribution of x 

Pp( X ) = ^[ 2 ^fi^da, (14) 
7T Jx V a - x z 

which can be computed numerically in the general case. 
By a linear change of variables T = 1— x/2 it provides the 
required distribution of transmission coefficient Pp(T). 
Moreover, for most interesting integer values of j3 the 
above integral can be evaluated analytically giving 

r ^K[y/T{2-T)] for (3 = 0, (CPE); 

PpAT)=l l\n 1+ ^- T) for (3=1, (COE); 

{ £T(2-T)D[y/T(2-T)] for (3 = 2, (CUE), 

(15) 
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where K and D stand for complete elliptic function of 
the first and the third kind, respectively H] . Somewhat 
more complicated derivation of this formula in the case 
of COE and CUE has already been given in where 
the elliptic function was expressed by the hypergeometric 
function as D(k) = §F(±, |;2;/c 2 ). 

We have generated numerically random unitary matri- 
ces of interpolating ensembles using the method proposed 
in |l(| . Figure [l] shows the distribution of transmission 
coefficients P(T) for M = 2 and the Poisson-COE transi- 
tion for three different values of the control parameter 5. 
Two narrow lines represent the formula ( |l5| ) with (3 — 
and (3—1, characteristic for the limiting cases of CPE 
and COE. The bold line in the figure b), obtained in the 
interpolating case, represents the best fit of (|l4| ) with 
(3 = 0.38 (a coincidence between the values of S and (3 
is accidental). Good quality of the fit reveals a certain 
validity of this formula with non-integer values of (3 for 
ensembles in between the usual universality classes. 




FIG. 1. Distribution of the transmission coefficient P(T) 
for Poisson - orthogonal transition for M = 2 obtained of 10 
random matrices. The control parameter S is equal to a) 0.0 
(CPE); b) 0.4; and c) 1.0 (COE). Solid and dashed narrow 
lines represent COE and CPE results, respectively, while the 
bold line in Fig. b) stands for the best fit with (3 = 0.38. 



C. Distribution Pp(T) for large M 

Since the transmission coefficient T is a function of M 
random variables with finite variances, one expects its 



distribution to be Gaussian for large M. This fact can 
be proved rigorously for CUE. For this case the traces 
z = tr(U) are distributed (in the limit M — > oo) as 
isotropic complex Gaussian variable |D| , what guaran- 
tees the Gaussian distribution as well for x =Re(z), as 
well for T — (M — x)/2. The Gaussian property also 
holds in Poissonian case, for which x is the sum of M 
independent terms, each being a cosine of uniformly dis- 
tributed random phases. 

For M >> 1 we conjecture the Gaussian distribution 
of P(T) for any value of (3 with the mean equal to M/2 
and the variance given by ( |To| ) . Numerical tests revealed 
Gaussian character of this distribution for any (3 already 
for M = 10. Figure || presents the distribution of trans- 
mission coefficient for the Poisson-CUE transition and 
three values of the control parameter 5. 
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FIG. 2. As in Fig. 1 for Poisson - unitary transition and 
M = 10. Solid lines represent the Gaussian distribution with 
(T) = M/2 and the variance given by @ for (3 = and 
(3 = 2 and found numerically for (3 = 0.4. 



V. CONCLUDING REMARKS 

Composed ensembles of random unitary matrices 
where applied for analysis of transmission coefficient in 
symmetric chaotic systems described by S-matrices with 
block symmetries. We simplified the derivation of distri- 
bution Pp{T) presented by Gopar et al. |)| and gener- 
alized their results proposing an family of interpolating 
distributions. 
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It should be noted that the Poissonian case j3 = dis- 
cussed above is not capable to describe an effect of local- 
ization. It corresponds rather to the case of scattering on 
a half transparent mirror, at which each scattering mode 
acquires a random phase shift. Another generalization of 
the model designed to describe effects of localization is a 
subject of a subsequent publication ]l~5| ]. 
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